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Abstract 

Robustness of biochemical systems has become one of the central ques- 
tions in systems biology although it is notoriously difficult to formally cap- 
ture its multifaceted nature. Maintenance of normal system function de- 
pends not only on the stoichiometry of the underlying interrelated compo- 
nents, but also on a multitude of kinetic parameters. Invariant flux ratios, 
obtained within flux coupling analysis, as well as invariant complex ratios, 
derived within chemical reaction network theory, can characterize robust 
properties of a system at steady state. However, the existing formalisms 
for the description of these invariants do not provide full characterization 
as they either only focus on the flux-centric or the concentration-centric 
view. Here we develop a novel mathematical framework which combines 
both views and thereby overcomes the limitations of the classical method- 
ologies. Our unified framework will be helpful in analyzing biologically 
important system properties. 
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1 Introduction 



Biochemical networks have evolved to operate in the face of internal and external 
perturbations (Kitano 2004). The response to these perturbations has shaped 



the systemic architectural blueprint comprising multiple layered and interrelated 
components (e.g., genes, proteins, metabolites). The dynamic processes involv- 
ing network-related biochemical components depend on a multitude of kinetic 
parameters, which remain elusive even for medium-size systems. Therefore, 
methods establishing a connection between structure and dynamics of biochem- 
ical systems hold the promise to enable the rigorous study of processes taking 
place on the underlying biochemical networks both at steady-state as well as 
dynamic setting. 

Two different classes of approaches have been developed to facilitate para- 
meter-independent analysis of biochemical networks: (i) flux-focused approaches 
including: flux balance analysis (FBA) (Varma and Palsson 1994) and its 
derivatives -flux variability analysis (FVA) (Mahadevan and Schilling 2003) 



and flux coupling analysis (FCA) ( Burgard et al. 2004 Marashi and Bockmayr 



2011 ), elementary flux modes (EFMs) (Schuster et al. 2000), and extreme path- 



ways ( Schilling et al. 1999); and {ii) concentration-centric approaches, rooted in 



chemical reaction network theory (CRNT) (Horn and Jackson 1972 Feinberg 



1979, 19951 and stoichiometric network analysis (Clarke 19881. 



Given a biochemical network, FBA relies on a linear programming formu- 
lation to calculate the steady-state fluxes under the assumption that the in- 
vestigated organism operates towards optimizing an objective function (e.g.. 



optimizing yield for metabolic networks (Varma and Palsson 1994)). FVA also 



has a linear programming formulation, with the aim of calculating the minimum 
and maximum values of individual steady-state fluxes for a particular value of 
the objective. FCA can be used to determine pairs of reactions whose flux ratio 
is the same in each steady state under the same environmental conditions. Like 
FBA and FVA, this approach can also be cast as a linear program. On the 
other hand, approaches based on EFMs allow decomposition of a given network 



into its smallest functional units operating in a steady state ( Schuster et al 



2000 Schilling et al. 19991. Although the problem of determining the set of 



all EFMs for a given biochemical network is computationally demanding, re- 
cent parallelized implementations of algorithms for EFM computation facilitate 



EFM-based analysis of genome-scale metabolic networks (Terzer and Stelling 



2008). Essential to both flux-based approaches is the usage of the underlying 



stoichiometric matrix which, without a specified kinetics, cannot be employed 
to make statements about steady-state metabolite concentrations. 

In contrast, CRNT uses mass-action formulation to study the qualitative 
behavior of the steady-state concentrations of the components regardless of the 
parameter values, i.e., for all steady-state reaction fluxes of the mass-action 
system satisfying the constraints imposed by the stoichiometry. The results 
of this framework answer questions related to the possibility for existence of 
multiple steady states, and rely on a structural index determined by interleaving 
the graph-theoretic and stoichiometric descriptions of the investigated network 
(Horn and Jackson 1972 Feinberg 1979 1995 Gunawardena^ ^2003 Conradi 



et al. 2007). 



Biochemical network invariants are of particular interest specifically because 
they relate to the principle of homeostasis. For instance, under the steady- 



state assumption, the concentrations of components do not change and, thus, 
are invariant. However, invariants in biochemical networks can be defined not 
only with respect to changes over time, but also changes with respect to different 
steady states that the system may assume under same environmental conditions 
(i.e., initial conditions and/or constraints). Note that the latter excludes the 
analysis of trivial invariants which are imposed in the form of conservation 
relations ( Schilling et al.[[T999[ Heinrich and Schuster 1996). 



In other words, invoking the steady-state assumption may induce additional 
invariants with respect to individual components or their combinations, which 
can ultimately reveal possible reduction in complexity of the system. As already 
stated, FCA provides the means for determining pairs of reactions whose ratio 
of fluxes is the same in each steady state the system may assume. 

In general, changes in fluxes and concentration, as key descriptors of the 
transitional behavior in biochemical networks, depend on each other. This 
stems from the fact that the reaction rate, i.e., flux, is cast as a function of 
the concentrations of the considered components. Therefore, the question arises 
whether there exist invariants on the level of concentrations and, if so, whether 
there is a connection between flux- and concentration-invariants. The answer 
to this question of course depends on the choice of kinetic law providing the 
relation between reaction fluxes and concentrations. 

Here, we focus on mass action kinetics, representing the simplest and most 
fundamental law of kinetics, to establish a connection between flux and con- 
centration-invariants. By interleaving the flux- and concentration-invariants, 
we provide a fundamentally new theoretical approach which can be used to 
uncover dependencies between fluxes and between concentrations, ultimately 
leading to a better understanding of system complexity. 

Therefore, our study establishes a connection between the two different views 
of computational systems biology — the flux-centric and the concentration-centric 
view. Since the theories and methods pertaining to the two views use different 
notations, a brief overview is provided to describe the used notation. 

2 Methods 

In chemistry, the law of mass action was established by Guldberg and Waage in 



the nineteenth century (Guldberg and Waage 1899 Abrash 1986). It assumes 



a mixture of large numbers of components which are homogeneously distributed, 
allowing approximation of the components' behavior with continuous variables. 
A reversible reaction, i.e., a reaction which can proceed in the forward and 
backward direction, is split into two reactions — the irreversible forward reac- 
tion and the irreversible backward reaction. The components consumed by an 
irreversible reaction are called substrates, while those produced are referred to 
as products. A reaction's rate is then modeled to be proportional to the product 
of the concentrations of the participating substrates, especially in the case of 
an elementary reaction which cannot be further divided into intermediate steps 
(Moore 1986 pg. 385). Under realistic chemical conditions, it is often the 
case that a given reaction almost certainly proceeds in one direction. In this 
situation, with the assumption that the reaction rate in one of the directions 
can be neglected, the reaction is treated as irreversible. Therefore, most mod- 
els of biochemical networks consist of a mixture of reversible and irreversible 



reactions. 

Here, for the application of specific theoretical methodology, each biochemi- 
cal network must be transformed to an equivalent one that consists only of irre- 



7?, U7?, 
y the subset 



versible reactions. Such a transformation is performed as follows (Gagneur and 
Klamt 2004 1 : Let the complete set of reactions be denoted by 7?. = 
where TZirr denotes the subset of irreversible reactions and TZr 
of reversible reactions. The set of reactions Ti-'^^^ is derived by splitting each 
reversible reaction from TZrev into two irreversible reactions, one in each direc- 
tion. The original network can then be described by a new set of reactions 
TZ' = TZ'^rr ^ T^irr with \TZ' \ — 2\TZrev\ + |^irr|- The Starting point for our 
methodologies derived here is always a biochemical network which is of this 
form, i.e., we assume that TZ denotes a set of irreversible reactions (see Example 



0- 



Example [TJ The eight irreversible reactions in the set 7?. — 

{Ri, R2, Rs, R4, R5, Re, R7, Rs}, given by 



1 


= A^X 


i?5 


= A+C^D 


2 


= X^ A 


Re 


D^A-hC 


3 


= B^X 


Ri 


= B + C^E 


4 


= X^B 


i?8 


E^B + C 



can in fact be regarded as four reversible reactions. The reversible re- 
actions are formed by Ri and R2, R3 and i?4, R^ and Rq as well as Rj 
and Ra,. 



The results from flux-centric approaches rely on investigating vector spaces 
associated to the stoichiometric matrix N (see Examplel2|. The principal object 
in the flux-centric approaches is given by the reactions and their fluxes. Here 
the term "flux" is used synonymously to "reaction rate" . A crucial vector space 
is that of the kernel of the stoichiometric matrix A'', which is represented by 
the set of flux vectors v fulfilling Nv = 0. Thus, the kernel of N describes all 
possible steady-state fluxes of the considered biochemical system. 



Example [2} The set of reactions from 
following stoichiometric matrix: 


Example 1 give rise to the 


^ CN 01 '^ lO 

D^ D^ C:^ D^ (i^ 


c? c^" c^" 






r -1 1 0-1 


1 0] 


A 




0-110 


-1 1 


B 


N = 


0-1 
1 


1 -1 1 
-10 


C 
D 







1 -1 


E 




1-1 1-1 





X 



The concentration-centric approaches, represented by CRNT, use a notation 



which combines linear algebra and set theory ( Gunawardena 2003 1 . For a given 



set of reactions, the set of complexes C is composed of the left- and right-hand 
sides of each reaction arrow. Any reaction y ^ y' G TZ can then easily be defined 
in terms of its complexes y, y' G C. 

Results from CRNT establish a relationship between the structure of a mass 
action system's network and the nature of the set of equilibria of the correspond- 
ing system of ODEs, independently of the rate constants ( Feinberg 1995 1 . Let 
P = {a; € M I a; > 0} be the set of positive real numbers. In the following, 
it is assumed that, if the system of ODEs of a mass action system admits an 
equilibrium, then the species' concentrations satisfy the following condition: 

Definition 111 Let c be the vector of concentrations of a mass action system. 
The system admits a positive steady state if dc/ dt = and c G P. 

A reaction network which mathematically captures the graph-theoretic prop- 



erties of a chemical mass action system is defined as follows (Feinberg 1979 



1995 Gunawardena 2003) 



Definition l2| A reaction network is a triple {S,C,TZ) where S is a finite 
set of species; C is a finite set of multisets of species, called complexes; TZ is 
a relation on C, denoted hy y -^ y' for y, y' G C, which represents a reaction 
converting y to y' . 



Based on a reaction network, a chemical reaction can be defined as follows 
( |Feinberg[ |1979[ |1995[ | Gunawardena] |2003| ) (see also Example [S]) : 



Definition Is} A ciiemical reaction network {S,C,TZ,IC) is a reaction 
network endowed with a function /C : 7?, — > P which associates a positive rate 
constant to each reaction of the reaction network. 



Example [3[ The set of reactions in Example [T] gives rise to the 
chemical reaction network {S,C,TZ,IC) with 


5 = {A,B,C,D 


X}, 


representing the set of species, 




C = {{A}, {B}, {A, C}, {B, D}, {D}, {E}, {X}}, 


representing the set of complexes, 




TZ ^= {i?i, R2, R3, Ri, R5 


Re,R7,Ra}, 


representing the set of reactions, and 




^ = {fcfli , fcfla : ^^3 ' f^Ri J ^fls 


■,kB.(iikR.T,kji^}, 


representing the set of rate constants. 





In the following, we use the concept of a partition of a set (see Definitions 
ql] - 93^ in the Appendix) to obtain new insights into the dynamics of chemical 
reaction networks. Of particular interest are the partitions of the set of reactions 
and those of the set of complexes (see also Example l4| , both of which can be 
placed in a relation defined as follows: 

Definition HI Given a set A and two equivalence relations ~ and ~' on the 
elements of A such that \/V' e A/ ~' 3V € A/ ^ with V CV. Then, A/ ~ 
is said to be a coarser partition than Aj ^' and Aj ~' is said to be a finer 
partition than Aj ~, denoted by Aj ~' < Aj ~. 



Example [4[ 

(a) The set 7^/ ~= {{i?i,i?2, ^3,-^4}, {^5,-^6,-^7,^8}} satisfies the 
conditions of Definition dsj Therefore, {i?i,i?2,-R3, -R4} and 
{-R5, Rq,Rt,Rs} are equivalence classes since i?i ~ i?2 ~ -R3 ~ -R4 
and i?5 ^^ Rq ^^ R-j ^^ Rg,. 

(b) The set 7^/ ~'= {{Ri,R2},{R3,R4},{R5,R6,R7,R8}} also sat- 
isfies the conditions of Definition ^ 

(c) The set {{i?i, i?2, -R3, -^4,-^5}, {-R5, -Re, -R/j-Rs}} does not satisfy 
condition (ii) of Definition 93] and, therefore, does not represent a 
partition of TZ. 

(d) By Definition |4] it is not difficult to see that 7^/ ~' < 7^/ ~. 



3 Results 

3.1 Invariant reaction ratios 

Consider two reactions y — >■ y', y — > y' G 7?. whose rate ratio at each positive 
steady state c is an invariant, i.e., [v{c)]y^y' /[v{c)]y^y' = const. Such pairs of 
reactions can readily be obtained by solving a linear program typical to FCA. We 
note that in FCA, a pair of reactions satisfying [v{c)]y^yi /[v{c)]y^y' = const is 
referred to as fully coupled reactions. Then, for two positive steady states c, c*, 
it must hold that [v{c)]y^y' /[v{c)]y^y> = [v{c*)]y^y' /[v{c*)]y-yy' whlch can be 
rewritten as [v{c)]y-fy' /[v{c*)]y^y> — [v{c)]y^y' /[v{c*)]y^y' . This leads to the 
following definition: 

Definition^ Let {S,C,TZ,IC) be a chemical reaction network that admits 
a positive steady state. Two reactions y ^ y',y ^ y' G 7?. are said to be in 
relation f.^ if for any pair of positive steady states c, c* , [v{c)]y^yi /[v{c*)]y^y' — 

[Vic)]y^y'/[vic*)]y^y>. 

The equivalence relation <^ induces a partition ofTZ into equivalence classes. 
The resulting partition is denoted by 7?./«^. For a given chemical reaction 
network, the partition 7?./<^ is an inherent property of the corresponding sys- 
tem of DDEs. The problem now is that of determining the partition 7?,/«^, 
i.e., all pairs of reactions whose flux ratio at a positive steady state is a con- 
stant. To address this nontrivial problem, we first observe that for two reactions 
y ~^ y'l y ~^ y' € 7?. with y = y, i.e., pairs of reactions which use the same sub- 
strate complex, it always holds that y — ^ y' <^y — > y', as shown by the following 
lemma: 

Lemma pi Let {S,C,TZ,IC) be a chemical reaction network. If two reactions 
y ^ y' ^y ^ y' ^ T^ share the same substrate complex, i.e., y = y , then 
y ^ y'^y-^ y' . 



Proof. The reactions y — > y', y — > y' G ^ share the same substrate complex. 
It follows directly that 

[Vic)]y^y> _ ky^y,Cy _ C^ _ fcy^S'^^ _ ^-'^^''^^V' 



Hc%^y, ky^y'c*y c*y ky^y,c*y Hc%^^,' 

and, therefore, y — ^ y' ^y ~^ y' ■ O 

Consider the partition of the set of reactions based on whether the reac- 
tions of the same equivalence class use the same substrate complex. Let the 
corresponding partition of the set of reactions be denoted by TZ/i^° . Lemma 
p^ shows that reactions, belonging to the same equivalence class in TZ/i^° , al- 
ways exhibit the same reaction rate ratio. Moreover, if two reactions use the 
same substrate complex, then they are element of the same equivalence class in 
7^/^. Then, for each 1?° G 11/^° there exists V G 11/^ with 2?° C V, and 
it follows that TZ/ ^^° < 7?./<^. The partitions of the set of reactions which are 
not finer than 7?./<^" and not coarser than 7?./<^ belong to the following set: 



Definitional Let {S,C,TZ,JC) be a chemical reaction network that admits a 
positive steady state. Then, 

Vn ■■= {7^/<^' I 7^/<^° < 7^/^' and 7^/<^' < 7^/^}. 

Clearly, the partition TZj J^ represents the finest element in Vfi- Moreover, 
it is also the finest partition of the set of reactions which is compatible with 
mass action kinetics, i.e., a partition such that all reactions with substrate 
complex y are in the same equivalence class. Analogously, the partition TZ/ <^ 
represents the coarsest element in Vn- If two reactions y ^ y\y ^f y' G 7?. are 
in the equivalence relation ^' with 7?./^' e Vui then it directly follows that 

[V{c)\y^y> /[V{c*)]y^yf = [v {c)]y^y- / [v {c* )]y^y- 'ATlA, thCrCfOrC, ?/ -^ J/'f^?/ "^ Tj' ■ 

The set Vn, together with the relation < on its elements, as given in Definition 
111 is a partially ordered set (see Lemma ^ in the Appendix for the proof) . 
Furthermore, together with the binary join and meet, from Definitions ^and 
^To|in the Appendix, V-r, represents a lattice (see Example [5] and Corollary ^14| 
in the Appendix for the proof). In contrast to the partition 7?./«^, the partition 
TZ/f^° can easily be determined by investigating the substrate complexes of the 
set of reactions. As a result, the higher a given element from V-r, is located in the 
lattice, the more information about invariant reaction ratios becomes available. 



Example [5[ The lower bound of the lattice corresponding to Example 
[3] is given by 

7^/^" = {{i?i},{i?2,i?4},{i?3},{i?5},{i?6},{i?7},{i?8}}, 

and, provided an oracle that yields TZ/i^, the upper bound is given by 

TZ/f^ — |{i?i, i?2, -R3: ^4}, {Rd^ReiRr, Rs}f- 

Then, V-ji consists of 71/^^° and TZ/f^ as well as all partitions of TZ 
that are coarser than TZ/i^" but finer than TZ/ 1^. 
(See the following examples for the derivation of this partition. For 
kn^ = 1, kii,^ = 1, kn^ = 1, ki^^^ = 1, fc_R^ = 1, k^^ = 1, fc_R^ = 1, 
kRg = 1, there exist at least two positive steady states: (i) ca = 1, 
cb = 1, cc = 1, CD = 1, ce == 1, cx = 1; and (ii) ca = 7/6, cb = 7/6, 
cc — 3/2, CD = 7/4, ce = 7/4, cx = 7/6; which implies that the 
rate ratio of Ri and Rg, is not constant for each positive steady state. 
Therefore, the partition consisting of two equivalence classes is in fact 
the coarsest.) 



One obvious question arises: What aspects of the chemical reaction network 
can be used to coarsen the finest partition in T'tj? It turns out that the kernel 
of the stoichiometric matrix can directly be used to answer this question, as 
illustrated in Example [6] (see also the concept of co-sets in Papin et al. (2004|). 



That is, the computation of fully coupled reactions from FCA using a linear 
program is not necessary. 



Example [6l A rational basis of the kernel of the stoichiometric matrix 
of Example [T] (see Example [2] for matrix N) is given by 



ker(iV) = 



1 








0] 


1 














1 











1 














1 











1 














1 











1 



Ri 

R2 

R3 

R4 

Rb 
Re 
R7 
Rh 



which implies the following partition of the set of reactions 
TZ/ ^— j^{i?i,-R2}, {Rs^RaJAR^^Rg}, {Rr.Rsjf- 

Since reactions R2 and R4 are in the relation ^^" , this implies that 
reactions i?i, R2, R3 and R4 satisfy Definitionis] Therefore, one obtains 

TZ/^" = -l^^{Ri, R2, R3, R4:},{R5i R6},{R7t Ra} f, 

with7^/«^" eViz- 



3.2 Invariant complex ratios 

Analogous to the previous subsection, consider two complexes y,y G C such 
that, at each positive steady state c, the ratio c^ j (9 is invariant. Then, for two 
positive steady states c,c* ^ it must hold that €}> jc^ = c*^/c*^, which can be 
rewritten as d> Ic*^ — 0^/0*^. This leads to the following definition: 

Definition^ Let {S,C,TZ,IC) be a chemical reaction network that admits a 
positive steady state. Two complexes y,y' £ C are said to be in relation «X> if 
for every pair of positive steady states c,c*, c^ j c*"^ — cP j c*"^ . 



The equivalence relation «^ induces a partition of C into equivalence classes. 
The resulting partition is denoted by C/<X>. Let the partition {{y} | y S C} of 
the set of complexes be denoted by C/^JU". Trivially, pairs of elements from 
an equivalence class in C / JL^° are, according to Definitional in the equivalence 
relation <~^. Furthermore, it is not difficult to see that C / ,^° < Cj ^S^. Con- 
sequently, the set of partitions of the set of complexes which are not finer than 
C/ ^° and not coarser than C/<X> can be defined, analogous to Vn^ as follows: 



Definition^ Let {S,C,TZ,IC) be a chemical reaction 
positive steady state. Then, 



ztwork that admits 



Vc :={C/^' I C/^" <C/^' andC/^' <C/^}. 



10 



The partition C/^S^ represents the coarsest element in Vc- Moreover, if 
two complexes y,y ^ C are in the equivalence relation ^', with C/^S^' S Vc, 
then it directly follows that cP jcP* = d> jc*^ and, therefore, y^S^y. The set 
Vc, together with the relation < on its elements is a partially ordered set (see 



Lemma £16 in the Appendix for the proof). Furthermore, together with the 



binary join and meet, as given in Definitions ^17| and ^21| in the Appendix, 



Vc represents a lattice (see Corollary £25 in the Appendix for the proof). In 



contrast to the partition C/<^, the partition C/,^" can be trivially determined. 
As a result, the higher a given element from Vc is located in the lattice, the 
more information about invariant complex ratios becomes available. 

3.3 Vtz is homomorphic to Vc 

In this section, it is shown that there exists a map ip : Viz ^^ Vc , which preserves 
the structure of V-ji- Let TZ/^^' G Vtz and V e TZ/^^' . Furthermore, let 
Uv — {y \ y ^> y' S 2?} and let U represent the set of complexes that do not 
participate in any reaction as a substrate. 



Definition 10, Let {S,C,TZ,IC) be a chemical reaction network that admits 
a positive steady state. Further, suppose that some partition TZj J^ G Vfi is 
given. The map ip : Vtz ^ Vc is defined by 

ifin/^' ):^{Uv\Ve n/^' }U{{y}\yeU}. 

In fact, one can show that Vtz is homomorphic to Vc with respect to ip (see 
Corollary £ 29 in the Appendix for the proof) . Example [7] illustrates the result 



of applying the map ip on the partition from Example [61 



Example [t} Partition 7^/^" = {{i?i, R2, R3, R4}, {-R5, Re}, 
{i?7, Rs}} from Example [6] can be mapped to a partition of the set of 
complexes: From Definition \T0\ it follows that 

ip(TZ/^") = \l^{Ri,R.2,R3,B4}T^{B.5,R6}T^{B.7,Bg}f, 

with U{B.i,R2,R.3,Ri} = ^A^X,X^A3^X,X^B} == {{A}, {B}, {X}}, 
^{R^,,Re} = W{A+C^D,D^A+C} = {{A, C}, {D}} and Z^{j^7,_Rs} = 

Z^{B+c^E.E^B+c} ={{B,C},{E}}. Therefore, 
^(7^/^")={{{A},{B},{X}},{{A,C},{D}},{{B,C},{E}}}. 



Now consider two positive steady states c, c* and a partition of the set of 
complexes Cj JL^ € Vc, for instance derived by (p from a partition TZ/^^' S 
Vtz- Since y,y e U with U £ C/^S^' satisfy 0^/0*^ = 0^/0*^, they also satisfy 
\a{cy/c*y) = ln(ci'/c*^) and, thus, {y - y)ln{c/c*) = 0. It is not diflicult to see 
that ln(c/c*) is orthogonal to any element in span({j/ — y \ y,y (z U}). 



11 



Lemma 11, Let (S,C,TZ,]C) be a chemical reaction network that admits a 
positive steady state. Furtherm,ore, consider a partition of the set of complexes 
C/<^' G Vc- If there exist y,y gC such that y — y£ span({j/' ~y' \ y' ,y' £ U}) 
with U e C/<X>' , then y^y. 

Proof. Let c, c* be two positive steady states of the chemical reaction network. 
The vector ln(c/c*) is orthogonal to each vector in span({y' — y' \ y',y' £ U}) 
and, thus, ln(c/c*) is also orthogonal to y~y. It follows that {y — y)\n{c/c*) = 
and, therefore, y^^S^y. □ 

Here we note that Lemma ITT] is very similar to Proposition S4.1 in |Shinar| 



and Feinberg (2010) (Supporting Information). Shinar and Feinberg use their 
proposition to investigate invariance of concentrations in chemical reaction net- 
works. However, with our lemma we attempt to coarsen the partition of the 
set of complexes. We show that the knowledge of the coarsest partition reduces 
the complexity of the system since dependencies between reactions and com- 
plexes become visible. In fact, by following the concentration-centered view of 



CRNT, as pursued in Shinar and Feinberg (2010), one may neglect the informa- 



tion contained in the dependence between reactions and complexes. Example [8] 
illustrates the implications of Lemma [TT| 



Example^ Lemma 11 applies to the partition C/<^" = tp{TZ/^") 
from Example [t] as follows: U = {{A}, {B}, {X}} e C/^" and y = 
{A, C}, y' - {B, C} e C. Then, y-y€ span({y' - y' \ y', y' e U}) since 
there exist y' , y' &U with y' -y' = {A}-{B} = {A, C}-{B, C} = y-y' . 
This results in the following, coarser partition of the complexes 

C/J^' = {{{A},{B},{X}},{{A,C},{D},{B,C},{E}}}. 



Let Vu :^ {y ^ y' \ y € hi and y —> y' G TZ} be the set of reactions which 
have the complexes in U as substrates. We next need the following definition 
which maps a partition of the set of complexes back to a partition of the set of 
reactions: 



Definition 30, Let (S,C,TZ,IC) be a chemical reaction network that admits a 



positive steady state. Furthermore, consider a partition C/^S^' £ Vc- The map 
fj, : Vc ^>^ Viz is defined by 

li{C/^'):={Vu \UeCI^'}\%. 

Application of the map /i to a partition of the set of complexes is illustrated 
in Example [9] 
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Example [gj Partition C/^' = {{{A}, {B}, {X}}, {{A, C}, {D}, 
{B, C}, {E}}} from Example [S] can be mapped to a partition of the set 
of reactions: From Definition [301 it follows that 



m(c/<. 



{^( 



{{A},{B},{X}},2?{{A,C},{D},{B,C}.{E}} 



}. 



with 2? 



{{A},{B},{X}} 



{i?l,i?2,i?3,-R4} and 2?{{A,C},{D},{B,C},{E}} 



{R5,Re,R7,Rs}- Therefore, 

^(C/^') = -{^{i?!, i?2, ^3,-^4}, {^5, ^6,^7,-^8}/ • 



It should be apparent that there may exist partitions in Vc that are coarser 
than (p{TZ/f^). Such partitions can consist of equivalence classes Z^ with \U\ > 2 
which also contain complexes from U. This information is lost when mapping 
such a partition of the set of complexes to a partition of the set of reactions 
using /i. As a result, ip is an injective map while /i is surjective. 

Consequently, the map ip has the useful property that a coarsening of a par- 
tition of the set of reactions implies a coarsening of the corresponding partition 
of the set of complexes. Furthermore, the map fi has the property that a coars- 
ening of the partition of the set of complexes guarantees that the corresponding 
partition of the set of reactions is not refined. 



3.4 Connections to CRNT 



Shinar and Feinberg (2010) defined an equivalence relation on the complexes, 
denoted by <~^, as follows: two complexes y^y d C are in the relation «^^ if, for 
two positive steady states c, c*, (y — y)ln(c/c*) = is satisfied. This equivalence 
relation induces a partition of the complexes. It is not difficult to see that this 
condition is equivalent to our Definitio n [gl 

Interestingly, Shinar and Feinberg (20101 specified which complexes are in 
the relation given by Definition |8] only for special classes of networks. CRNT 
defines several properties of reaction networks, e.g., the deficiency and weakly 
reversibility (Feinberg 1979 1995 Gunawardena 2003). The deficiency 6 is 



an index which can provide information about the dynamic behavior of a mass 
action system independently of the rate constants. This index is computed 
as (5 = n — I — q, where n is the number of complexes, I is the number of 
linkage classes, i.e., the number of connected components of the graph which 



can be build by the reactions (see Example 10 for an illustration), and q is the 
rank of the stoichiometric matrix TV. A reaction network is said to be weakly 
reversible, if, in each connected component, there exists a path from any node 
(i.e., complex) to all other nodes in the connected component. For instance, 
in deficiency-zero reaction networks which are weakly reversible, complexes are 
in the same equivalence class if they are elements of the same linkage class 
(Shinar and Feinberg 2010). Furthermore, in deficiency-one reaction networks. 



all complexes of nonterminal strong linkage classes are elements of the same 



equivalence class (Shinar and Feinberg 2010). With our definition of the lattice 



Vc, it is clear that in both cases the resulting partitions are elements of Vc- 
Nevertheless, they may still not represent the coarsest element. 
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In contrast to the results from CRNT, our findings do not pertain to special 
classes of networks. In fact, our theoretical results deal with equivalence classes 
of reactions which can be mapped to a partition of the set of complexes in 
general networks. The defined map in turn enables the determination of coarser 
partitions of the set of complexes from coarser partitions of the set of reactions. 



Example |10[ The chemical reaction network from Example [3] can be 
represented by the following graph: 

A T^ ^ X ^ ^ R 

A+C i — ^ D 

kRg 

feny 

B+C i^ ^ E 



The set of nodes of this graph is equivalent to the set of complexes (here 
n — 7). The connected components of this graph are equivalent to the 
linkage classes (here / = 3). The strongly connected components of this 
graph are equivalent to the strong linkage classes (here equivalent to 
the linkage classes). The terminal strong linkage classes are the strong 
linkage classes from which there exists no path to another strong linkage 
class (here equivalent to the linkage classes). Thus, this reaction network 
is weakly reversible. The rank of the stoichiometric matrix N is q = 
4 (see also Example [2]). As a result, the deficiency of this reaction 
network is 6 — n — I — q — 0. Consequently, the coarsest partition of 
the set of complexes with respect to Definition [8] contains at most three 
equivalence classes. However, by applying our approach, there exists a 
coarser partition consisting of two equivalence classes, as illustrated in 
Example |9] 



4 Conclusions 

We analyzed ratio invariants in chemical reaction networks on the level of re- 
actions and on the level of complexes. We show that there can exist pairs of 
distinct reactions whose reaction rate ratio is constant in each steady state. 
This fact can be used to define a partition of the set of reactions. The knowl- 
edge of all such pairs of reactions determines the coarsest partition of the set 
of reactions (7?./<^). Furthermore, we defined the set of partitions which are 
compatible with mass action kinetics {V-r), i.e., reactions which use the same 
substrate complex must be elements of the same equivalence class. As a result, 
the coarsest partition of the set of reactions cannot be refined arbitrarily with- 
out violating constraints imposed by mass action kinetics. The idea of existence 
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of the coarsest and finest partition {TZ/^^") of the set of reactions, led to the 
introduction of the lattice of partitions of the set of reactions, represented by 

Vn- 

Analogously, we defined a finest (C/^^") and a coarsest partition (C/«^) 
on the set of complexes as well as the corresponding lattice Vc ■ Since the rate 
of a given reaction is directly influenced by the product of the concentrations 
of its substrates, we introduced an injective map (tp) which converts a partition 
of the set of reactions to a partition of the set of complexes. Additionally, we 
defined a surjective map (/x) which converts a partition of the set of complexes 
to a partition of the set of reactions. 

We explicitly point out that the coarsest partition of the set of reactions, 
and, thus, of the complexes, might depend on the rate constants. As illustrated 
by the examples, there always exist partitions of the set of reactions which can 
be derived independently of the rate constants. However, at this point, it is not 
clear whether the coarsest partition of the set of reactions changes for different 
choices of the set of rate constants. 

The study of metabolic networks has been hampered by the dichotomy in 
the computational approaches focused either on analysis of biochemical network 
fluxes or on concentrations of biochemical network components. Our study 
bridges this dichotomy by establishing the relation between the flux-centric and 
concentration-centric approaches, here represented by FCA and CRNT. Based 
on the established connection between invariants on the level of reactions and 
on the level of concentrations, we provide a method that will allow a deeper 
algebraic insight into the dynamic behavior of chemical mass action systems 
avoiding numerical computations. Our theoretical findings also provide the 
impetus for rigorous analysis of biological systems concerning the dynamics of 
biochemical networks. 
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A Appendix 

A.l Prerequisites 

The following general definitions are indispensable to understand the notation 
we use in this work: 

Definition ^l] A relation ^ on a set A is called an equivalence relation, 

if it satisfies 

(i) X '^ X \ Vx G A (reflexivity), 

(ii) if x,y Cz A and x ^ y, then y ^ x (symmetry), and 
(Hi) if x,y, z £ A, X ^ y and y ^ z, then x ^ z (transitivity). 

Definition ^2} Consider a relation ^ on a set A and an element x E A. The 
set of elements y £ A with x ~ y is called the equivalence class of x. The set 
of equivalence classes of A is denoted by A/ ~. 

Definition ^3} Let A be a set and letVi,. . . ,2?„ C A. The set {T>i, . . . ,2?n} 
is called a partition of A if and only if 

(i) ViU ...UVn^ A and 

(ii) V, n Pj = with V,,Vj £ A and Vi ^ Vj . 

The equivalence classes of a set A with respect to a relation ~ yield a parti- 



tion of ^ (Makinson 2008). Therefore, we use the terms "equivalence classes" 



and "partition" of a set with respect to a given relation synonymously. 

Given a chemical reaction network {S,C, TZ, JC) and a partition Ti./ ~ of the 
set of reactions, we define U-d := {y \ y ^ y' ^V} for V G TZ/ ~. Furthermore, 
let _ 

U:^C\ IJ Uv 

denote the set of complexes that do not participate in any reaction as substrates. 
Additionally, let y := {{y} \ y G U}. Analogously, given a partition C/ ~ of 
the corresponding set of complexes, we define Du := {y ^>- y' \ y £ U and y — >■ 
y' G TZ} for U €C/ ^. 



•R 
4- — -^ 



For graph-theoretical concepts, see Bollobas (1998). 

A. 2 Viz is a lattice 

The set T'-jz is defined as follows (see the main text for the definition of TZ/ 
andTZ/^): 

Definition 94} Let {S,C,TZ,IC) be a chemical reaction network that admits a 
positive steady state. Then, 

Vn := {TZ/^' \ TZ/^" < 7^/^' and TZ/^' < TZ/^}. 
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Then, Vn is a partially ordered set as the following lemma shows: 

Lemma 95| Let {S,C,TZ,IC) be a chemical reaction network that admits a 
positive steady state. The set Viz is a partially ordered set. 

Proof. For A, A' , A" G Pk, it holds that: 

• A < Ahy definition. 

• liA<A! and A! < A, then \fV € A, 3V' e A' with VCV' and WV £ A', 
3V" e A with V C V" which implies VCV". Since V, V" e A, it follows 
that V = V", implying V QV' &nAV' QV. Therefore, A = A'. 

• liA<A' and A' < A", then VP e A, 3V' e A' with PCX*' and VP' G 
yl' 32?" e yl" with V C X>". It follows directly that yV e A 3V" e A" 
with 2? C V" and, therefore, yl < ^". 

Therefore, Vn is a partially ordered set. D 

The following argument will be crucial in what follows: Let A,A\A" be 
partitions of TZ and let A < A" and A' < A". Furthermore, let V £ A and 
V G A'. From the properties of a partial order (see Lemma ^5]), w.l.o.g., it 
follows that 3V'{,V'^ e yl" with V C 2?'/ and V C V'^ lfVnV'^9, then 
2?'/ = 2?'2', since A" is a partition of 2^, implying VUV CV'{. 

Definition am Let {S,C,TZ,IC) be a chemical reaction network that admits a 
positive steady state. The binary operation AAA' of two elements A, A' € V-jz 
IS defined by A A A' := {V D V \V gA and V £ A' and VnV ^ 0}. 



Lemma S[7l Let {S,C,TZ,IC) be a chemical reaction network that admits a 
positive steady state. Given two elements A, A' G T'k, the set AAA' is a 
partition of the set of reactions. 

Proof. For each y —>■ y' G TZ, there exist exactly one V G A and exactly one 
V £ A' with y —i' y' £ V and y —>■ y' £ V , since A and A' are partitions of the 
set of reactions, so that y —^ y' £ V \^V' ^ % . It follows directly that for each 
y ^ y' £ TZ there exists exactly one V£AaA' with y ^ y' £ V. As a result, 
the union of all elements in AAA' equals 2^ and 2? n 2?' = for V,V' £AaA' 
and V y^ V . Therefore, .4 A .4' is a partition of the set of reactions. D 



Lemma aSl Let {S,C,TZ,IC) be a chemical reaction network that admits a 
positive steady state. Given two elements A, A' £ Vu, the set AAA! is an 
element ofViz- 

Proof, liy ^ y' £V with V £ A, then Vf^y} C V. It follows that for each 
T^{y} ^ ^ there exists exactly one V £ A' with V^y^ CVO V . As a result, for 
each V^yj with y £ C \U there exists exactly one V" £AaA' with 2?{j,} C V" . 
Then, AAA' is at least as coarse as {2?{j,} \ y £ C\ U} which is in turn the 
finest element in V-r,- 
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Lemma ^ shows that AAA' is a partition of TZ. From Definition ^ it 
follows that yV" GAAA',3VeA with V" CV. Analogously, VP" eAAA', 
3V' e A' with V" C 2?'. Then, it also follows A A A' < A and A A A' < A' . 
Therefore, AAA' eVn- □ 

Lemma 99} Lei {S,C,TZ,IC) be a chemical reaction network that admits a 
positive steady state. If A, A' G Viz, then AAA' is the greatest element in V-r 
which is finer than or equal to A and A! . Therefore, the binary operation AAA' 
is the meet of A,A'. 

Proof. The statement is proved by contradiction: Suppose that 3A" with 
AAA' <A", A" < A and A" < A'. Further suppose that AAA' ^ A". Then, 
3V e A,V e A' and V" e A" with ij) y^ V nV G A A A' and V nV C V". 
From A" < A and A" < A' it also follows that V" C V and V" C V. But this 
directly implies V" CVn V which is a contradiction to P n P' C P". D 



Definition ^10, Let (S,C,TZ,IC) be a chemical reaction network that admits 
a positive steady state and let A, A' € Vn- Furthermore, let G := (V,£) be 
a graph with V = {V U V \ V e A, V e A', and V nV ^ $} and £ ^ 
{{v,v') \ v,v' eV andvTw' ^0}. LetG be decomposed as G = G'^^^ U .. .UG^"^ 
where G^*' = {V^^',£^^') with 1 < i < s represent the connected components. 
Then, the binary operation A\/ A' is defined by 

Aw A' := < [j w|l<i<sl. 



Lemma S |ll[ Let {S,C,TZ,IC) be a chemical reaction network that admits 
a positive steady state. Given two elements A, A' G T'k, the set AW A' is a 
partition of the set of reactions. 

Proof. Since two nodes of G are adjacent if the intersection of the corre- 
sponding sets is nonempty, it follows that y — >■ y' G 7?. is element of exactly one 
V G Ay A'. As a result, the union of all elements in Ay A' equals TZ, VnV = 
for V,V' e Ay A', and X> 7^ P'. Therefore, Ay A' is a partition of the set of 
reactions. D 



Lemma S12, Let (S,C,TZ,IC) be a chemical reaction network that admits a 
positive steady state. Given two elements A, A' G V-jz, the set Ay A' is an 
element ofV-jz- 



Proof. Lemma £11 shows that Ay A' is a partition of TZ. From Definition 
^0\ it follows that WV G A, 3V" G ^ V ^' with V C V" . Analogously, the 
same holds true for the elements of A'. Then, it follows that A < Ay A' and 
A' < Ay A' such that Ay A' is as coarse as the finest element in T^-ji. 

Let A" be the coarsest element in T'-ji, i.e., A < A" and A' < A" . If 
VnV ^H) with X» G ^ and P' G ^', then 3V" G A" with V C V" and V C V" 
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such that V\JV' CV". Now let ULi (^* U I?-) be an element oi AV A' . Then, 
3V'l, . . . ,I?i' e yi" with Vi U P'l C V'l, ...,VtVJV[C V'l. Let i, j e {1, . . . , t\ 
with i ^ j. liViUV^ and VjUVj are connected by an edge in G, then I?" — V- 
because A" is a partition of 7?,. Then, since an element oi A\/ A' represents 
a connected component of G, it follows that V'l = . . . = I?". As a result, for 
each U*=i(^i U I?9 e ^ V y^' 32?" G yl" with ULi(^« ^ 2?^) C 2?" and, thus. 
Ay A' < A". Therefore, AV A' eVn- □ 



Lemma 513, Let {S,C, TZ,)C) be a chemical reaction network that admits a 
positive steady state. If A, A' G Vn, then ,4 V ,4' is the smallest element in V-r. 
which is coarser than or equal to A and A! . Therefore, the binary operation 
Ay A' is the join of A,A' . 

Proof. The statement is proved by contradiction: Suppose that 3A" with 
A < A" , A' < A" and A" < Ay A'. Further suppose that Ay A' ^ A". 
W.l.o.g., this directly implies that 3 lJ*^^(2?jU2?^) G Ay A' with Vi, ... ,'Dt G A 
and V[, . . . ,V[ e A' a,nd t > 2 such that 3V" G A" with 2?" C ULi(^« ^ 2?^. 
Then, 3i G {1, . . . , i} with 2?, U 2?^ C 2?". Because ULi(^* U V[) is the union 
of all nodes of a connected component of G, 3j G {1, . . . , t} and i y^ j such that 
{Vi U 2?^) n (2?j U 2?^) 7^ 0. It follows that 2?, U 2?^ U 2?^ U V'j C 2?" since the 
sets represent equivalence classes of partitions. Repeating this argument for all 
nodes of the corresponding connected component results in [Ji^i{ViUV^) C V" 
which is a contradiction to V" C Ui=i(^j '-^^D ^i^"^; thus, to the assumption. D 



Corollary £14, Let {S,C,TZ,IC) be a chemical reaction network that admits 



a positive steady state. The set V-r is a lattice. 

Proof. Lemma ^ shows that Vn is a partially ordered set. Lemmas ^ ^ 
and ^ show that the set defined in Definition ^ is the meet of two elements 
in Viz- Furthermore, Lemmas ^TT] ^l2| and ^T3| show that the set defined in 



Definition £ 10 is the join of two elements in Vu- Therefore, V-r is a lattice. D 



A. 3 Vc is a lattice 

The set Vc is defined as follows (see the main text for the definition of C/ 



C 



and C/< 



c 



Definition £15 Let (S,C,TZ,IC) be a chemical reaction network that admits 



a positive steady state. Then, 

Vc :={C/^' I C/^° <C/^' andC/^' <C/^}. 



Lemma £16. Let {S,C,TZ,IC) be a chemical reaction network that admits a 



positive steady state. The set Vc is a partially ordered set. 

Proof. The proof works analogously to the proof of Lemma ^ D 
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Definition ^17[ Let {S,C,TZ,JC) be a chemical reaction network that admits 
a positive steady state. The binary operation AAA' of two elements A, A' S Vc 
IS defined by A AA' -.^ {U nW \U eA and W £ A' and UnW y^Hl}. 



Lemma S18, Let {S,C,TZ,IC) be a chemical reaction network that admits 
a positive steady state. Given two elements A, A' G Vc, the set AAA' is a 
partition of the set of complexes. 

Proof. The proof works analogously to the proof of Lemma 371 D 



Lemma S19, Let {S,C,TZ,]C) be a chemical reaction network that admits 
a positive steady state. Given two elements A, A' G Vc, the set AAA' is an 
element of Vc • 



Proof. Lemma E 18 shows that ^ A ^' is a partition of C. The partition 
{{y} \ y ^ C} is the finest partition of C and is element of Vc, i.e., {{y} \ y G 
C} < A A A' . The second part of the proof works analogously to the second 
part of the proof of Lemma ^ D 



Lemma ^20, Let {S,C,TZ,IC) be a chemical reaction network that admits a 
positive steady state. If A, A' G Vc, then AAA' is the greatest element in Vc 
which is finer than or equal to A and A' . Therefore, the binary operation AAA' 
is the meet of A, A'. 

Proof. The proof works analogously to the proof of Lemma ^ D 



Definition ^21, Let {S,C, TZ,]C) be a chemical reaction network that admits 
a positive steady state and let A, A' S Vc. Furthermore, let G := {V,£) be 
a graph with V = {U U W \ U <E A, W e A' , and U W ^ ?)} and £ = 
{{v,v') \ v,v' eV andvTw' ^0}. LetG be decomposed as G ^ G^^' U. . .UG^"*) 
where G'*' — (V'-*\f '*■') with \ < i < s represent the connected components. 
Then, the binary operation AW A' is defined by 

Ay A' :={ I I v\\<i<s\. 




Lemma S22, Let {S,C,TZ,)C) be a chemical reaction network that admits 
a positive steady state. Given two elements A, A' € Viz, the set A\/ A' is a 
partition of the set of reactions. 



Proof. The proof works analogously to the proof of Lemma £ 11 D 
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Lemma 5 23, Let {S,C,TZ,]C) be a chemical reaction network that admits a 
positive steady state. Given two elements A, A' G T'tj; the set A \/ A' is an 
element of Vc ■ 

Proof. The proof works analogously to the proof of Lemma 312] D 



Lemma 5 24, Let {S,C,TZ,IC) be a chemical reaction network that admits a 
positive steady state. If A, A' G Vc, then A\' A' is the smallest element in Vc 
which is coarser than or equal to A and A! . Therefore, the binary operation 
A\J A! is the join of A, A'. 



Proof. The proof works analogously to the proof of Lemma £ 13 D 



Corollary S25, Let {S,C, TZ,]C) be a chemical reaction network that admits 



a positive steady state. The set Vc is a lattice. 



Proof. The proof works analogously to the proof of Lemma £ 14 D 



A. 4 Viz is homomorphic to Vc 

In this section, we show that there exists a map (p : V-jz — > Vc that preserves 
the structure of Viz . 



Definition E26, Let (S,C,TZ,IC) be a chemical reaction network that admits 
a positive steady state. Further, suppose that some partition A G Vtz is given. 
The map if : Vu — > Vc is defined by 

ifiiA) := {Uv\VeA}U {{y} \ y eU}. 

First, it is shown that (p{AAA') = ip{A)Aif{A') for two elements A, A' G Vn- 
Let V e A, V e A' and A, A' e Vn- It is not difficult to see that U-vn-D' = 
{y \ y ^ y' e V nV'} ^ {y \ y ^ y' e V} n {y \ y ^ y' € V'} = Uv nUv 
Since y— >y'G7?.,j/— >y'Gl? implies y — > y' G P, it follows that V = T>u^. 
Consequently, X> n X>' 7^ if and only if Uvr\Uv' ^%. 

Lemma 9271 Let {S,C,TZ,IC) be a chemical reaction network that admits a 
positive steady state. Then, for two elements A, A' G Viz, it holds that ip{A A 
A') = ip{A)Aip{A'). 

Proof. From the previous arguments, it follows that 

(fiAAA') ^{Uvnv I VeA,V' £A',VnV' y^<ii}uy _ 

^{UvnUv \veA,v' eA',vnv' j^iDyuy _ 

= {UvnUv I Uv G ip{A),Uv' G ip{A'),Uvriiiv' ^9}uy 
^{A)Af{A') = {ur\U' \u e^{A),u' eip{A'),ur\U' ^9}, 

which proves the lemma. D 
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Second, it is shown that (p{A V A') = fiA) V fiA') for two elements A, A' £ 
Vn- It is also not difficult to see that U-vuv — {y \ y ^ y' £ ViiV} = {y \ y ^ 
y' €V}U{y\ y^y' €V'}=UvUUv 



Lemma S |28[ Let (S,C,TZ,IC) be a chemical reaction network that admits a 
positive steady state. Then, for two elements A, A' £ V-r,, H holds that ip{Ay 

A')^^{A)yip{A'). 

Proof. Let G-ji :— {VTz,£n) be the graph analogous to the graph in Definition 



£10 consisting of s-n connected components. Further, let Gc ■— {Vc,£c) be 
the graph analogous to the graph in Definition ^21[ consisting of sq connected 
components. Then, 

<f{A V A') - {Wu„^^(.) v\l<^<s^^}^y 

= {Upup"ev<,') ^^uP' \l<t<sn}uy ^ 
= iUpuD'eviJ' l^v^Kv \l<i<sn}uy 

ifiiA) V ^{A') = {U,„ev<*' « I 1 < * < 5c}, 
which proves the lemma. D 



Corollary £29, Let {S,C, TZ,]C) be a chemical reaction network that admits 



a positive steady state. The map Lp : V-r -^ Vc is a homomorphism. 



Proof. Lemmas £27 and £28 show that '.p{A A A!) = ip{A) A ip{A!) and 



(p{A y A') = ^{A) V ^{A'). Therefore, (p preserves the structure of V-ji- D 

Analogous to map ^p, there also exists a map n : Vc ^ Vu which converts 
partitions of the set of complexes to partitions of the set of reactions. 



Definition 30, Let {S,C,TZ,IC) be a chemical reaction network that admits 
a positive steady state. Furthermore, suppose a partition A £ Vc- The map 
jjL :Vc ^ V-ji, is defined by 

tx{A) :- {Vu\U£A}\ 0. 

While the map (/? preserves the structure of Vu: this is not true for /z. In 
the following we show that /i is a surjective function. 



Lemma S |31[ Let {S,C,TZ,K.) be a chemical reaction network that admits a 
positive steady state. Then, for two elements A, A' £Vc, it holds that ^{A /\ 

A') = ^i{A)^^l(A'). 

Proof. There are three different cases for Vunw and T>u n T>u' , respectively, 
where Z^ £ A.W £ A': 



22 



{\)uc^u' cu^ Vunw =Vun Vw = 0. 

(2) {UnU')nU ^9 (and not case (1)) => Vunw ^T^u^ T^w + 0- 

(3) {U^U')r\U = %^ Vunw ='Dur\ Vw ^ 0. 

That is, for cases (2) and (3), Vu G iJ-{A) and Vw G ^-i{A'). Then, 

Ijl{AaA') = {Vunu' \U eA,U' eA'.UnW ^9}\9 

^{VunVw \U eA,U' eA\unu' ^i!}}\i!} 
fi{A) A ^iiA') ^{vnV \ve ^iiA),v' e fiiA'),v nV ^ 0}, 

which proves the lemma. D 

For the remaining proof we need some additional definitions. Let Vc be the 
set of nodes in Gc with Vc ■— {v eVc | w C U}. Analogously, let £c be the set of 
edges in Gc with f c := { (w, v) e £c \ uHv C U}. Finally, let Gc := {Vc\Vc,£c\ 
£c) = {VcJc) = {&c\...,G'^c^} with G^*^ := (V^*\4''') being the connected 
components of Gc ■ It is not difficult to see that Gc represents a partition of the 
set of substrate complexes. Since each substrate complex corresponds to a set of 
reactions which are elements of the same equivalence class of a partition in V-jZj 
such a partition of the set of substrate complexes is equivalent to a partition 
of the set of reactions. Then, it directly follows that TZ/^^" < TZ/^^' , where 
7^/^' = M({U,ev('> I 1 < ^ < sc}) and TZ/^" = Ai({U„ev(o I 1 < * < t}). 



Lemma S32, Let {S,C,TZ,)C) be a chemical reaction network that admits 
a positive steady state. Further, let two elements A, A' S Vc be given such 
that, based on A\l A! , the graph Gc is defined and, based on IJ.{A) V fi{A'), the 
graph G-jz is defined. If TZj ^' = MdU^cy*'' I 1 ^ * ^ 0) ^'^'^ TZj ,^" — 
{U„pv'" I 1 < * < «k}, then 7^/^" = 7^/^"' . 

Proof. The proof consists of two parts: 

(1) ^ 

(l.i) U\JU' (^Vf ^UC\U' i-%%U^Vur\Vw -h^^^i with Vu u 
Vw e V^^ since U ^ A,U' (^ A! implies Vu E fi{A),Vw E niA'). 

(i.ii) {{u u W), {u u u')) e 4'^ ^ (u u W) n (u u u') 7^ g, u ^ 

(Vu UVw)n [V^ \JV^,)^%^ 3j with {{Vu UVw), {v^ u %,)) e 

(l.iii) Finally, from (l.i) and (I.ii) follows {Vuuw \UUU' E V^'^} C V^'^ 

(2) ^ 

(2.i) PUP'e V^^ ^VnV ^(d^UvnUv ^ $ ^ 3k with UvUUv e 
V^''^ since V e //(y^),!?' G /-i(yt') implies Uv e ^,Z^d' e A'. 

(2.ii) ((x>ux>'),(^u2?')) e^K^ ^ {vuv')n{vuv') 7^0^ {UvUUv')n 

{Uf,UUf,,) ^ ^ %U ^ 3k with ((Z^p U Z^x)')> (^D U Uf,,)) e 4'"'* =^ 
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;(j) <- fn \ 1 1 \ \i II r- -\^^^■)^ 



(2.iii) Finally, from (2.1) and (2.11) follows V^^ C {Vu^jw \U\JU' e V}^'}. 
Then, since Gq Is equivalent to a partition of the set of reactions, from (I.IU) and 



(2.111) It follows that V^'^ = V^''\ Therefore, Vi e {1, ... ,0 3j e {1, . . . , stj} 



(i)x _^,U} 



(and vice versa) with {Vu^w \UUU' eV^'}^ V!j(' . D 

Corollary ^33[ Let {S,C,TZ,IC) be a chemical reaction network that admits a 
positive steady state. Then, for two elements A, A' £ Vc, it holds that ^{A) V 

^x{A')<^l{AyA'). 
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shows that ^i{A)\/ ^i{A') = m({U pv'-' I 1 ^ « ^ 0)- As a 



c 



Proof. Lemma £ 

result, ^l{A) V //(ir- /^({U„6v<-' I 1 < ^ < 0) < /^({U„ev(*' I 1 < ^ < ^c}) = 

^l{Ay A!). " " u 
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